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I.  INTRODUCTION 


Of  all  the  elements  required  for  a  numerical  simulation  of 
seismic  coupling  from  a  nuclear  explosion,  the  most  important  and 
the  most  uncertain  are  the  physical  (constitutive)  models  of  each 
nonlinear  process  induced  in  the  medium  by  the  propagating  stress 
field,  and  the  material  properties  required  by  the  constitutive 
model  which  makes  the  simulation  site  specific.  These  two  elements 
cannot  be  validated  by  numerical  experiments.  In  fact,  there  is  a 
trade-off  between  the  two,  with  an  incorrect  constitutive  model 
requiring  incorrect  material  properties  in  order  to  produce  results 
which  agree  with  measured  ground  motion  data. 

Our  approach  to  resolving  this  uncertainty  has  been  to  base 
the  development  and  formulation  of  each  constitutive  model  on  the 
results  of  quasi-static  laboratory  tests  on  rock  samples,  and  then 
to  apply  the  model  to  the  simulation  of  explosion  induced  ground 
motion.  If  the  application  appeared  promising  in  terms  of  its 
ability  to  match  ground  motion  data  with  the  constitutive  model 
using  independently  measured  material  properties,  then  the  model  was 
retained  for  use  in  future  ground  motion  simulations. 

This  approach  has  shown  that  two  specific  constitutive  models 
severely  affect  seismic  coupling.  They  are,  the  irreversible 
collapse  of  air-filled  porosity  and  the  modification  of  material 
strength  due  to  pore  fluid  pressure,  i.e.,  effective  stress. 
Therefore,  the  validation  and  testing  of  these  models  are  issues  of 
continuing  importance  and  concern. 

In  this  report,  we  present  a  rather  severe  test  of  these 
models  by  using  them  to  simulate  experiments  ir  which  small  scale 
explosions  were  detonated  in  grout  spheres.  High  quality, 
reproducible,  particle  velocity  data  were  obtained  from  these 
experiments,  de  show  that  an  effective  stress  law  coupled  with  the 
irreversible  collapse  of  air-filled  porosity  provide  a  very  simple, 
straightforward  explanation  of  this  data.  These  results  have  served 
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II.  CONSTITUTIVE  MODELS  FOR  POROSITY  AND  EFFECTIVE  STRESS 


2.1  POROSITY 

The  first  laboratory  tests  on  the  irreversible  collapse  of 
air-filled  porosity  in  rock  samples  were  reported  by  Stephens  and 
Lilley  (1970)  in  which  they  summarized  their  work  on  rocks  of 
interest  to  the  Plowshare  Program.  Based  on  this  work,  a 
constitutive  model  for  porosity  collapse  was  formulated  (Cherry  and 
Peterson,  1970),  refined  (Cherry,  et_  aL,  1973),  and  used  to  match 
ground  motion  data  from  explosive  sources  (Cherry  and  Peterson, 
1970;  Riney,  et_ak ,  1973). 

These  results  were  encouraging  enough  to  suggest  that  we 
should  use  the  model  to  determine  the  effect  of  air-filled  porosity 
on  teleseismic  magnitudes.  Cherry,  Rimer,  and  Wray  (1975)  performed 
the  parameter  study  and  found  that  the  introduction  of  porosity 
causes  a  severe  reduction  in  magnitude  with  magnitude  changing  by 
0.17  units  when  porosity  changes  from  zero  to  one  percent. 

The  porosity  model  developed  by  Cherry,  et  (1973)  has 
essentially  remained  unchanged  for  our  applications.  In  the  model, 
the  hydrodynamic  component  of  the  stress  tensor,  P,  is  dependent  on 
air-filled  porosity  via  a  parameter,®,  given  by 


a 


+  V 


+  V 


u 


_v 

A 

V 


(2.1) 


where  Vr,  V^,  and  Vp  are  the  respective  rock,  water  and  air 
pore  volumes,  and  v  is  the  specific  volume  of  the  partially 
saturated  mix. 

The  material  properties  required  by  the  model  are: 

A  A 

1.  P(v),  the  pressure/volume  relation  for  the  rock  with 
air-voids  removed. 

2.  P  ,  the  pressure  at  which  all  air-voids  are  removed. 

3.  Pg,  the  pressure  above  which  irreversible  removal  of 
air-voids  begins. 
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4.  4>  ,  the  initial  air-void  volume  fraction,  where 

o  o 

and  a  are  related  by 
o 


a 


0 


(2.2) 


5.  kQ,  the  zero  pressure  bulk  modulus  of  the  porous 
mixture. 

Given  the  above  material  properties,  then 


P  -  i  P(v/«) 


(2.3) 


where 


»0  >  a  >  1  (2.4) 

and  the  material  properties  (P  .  P  .  <J>  ,  and  k  )  determine 

c  e  o  o 

the  variation  of  a  with  v  under  the  constraints  for  pore  collapse 
assumed  by  Cherry,  £t  aK  (1973). 

As  noted  earlier,  this  model  has  be^n  used  successfully  in  a 
variety  of  ground  motion  simulations  since  its  development.  In  this 
report  we  present  additional  evidence  that  the  model  adequately 
matches  observed  ground  motion  in  rock  materials  whose  air-void 
porosities  range  between  1.5  percent  and  13.4  percent. 

2.2  EFFECTIVE  STRESS 

The  development  of  constitutive  models  for  an  adequate 
treatment  of  material  strength  has  been  an  issue  of  debate  and 
uncertainty  among  the  ground  shock  calculation  community  for  over 
ten  years.  The  major  issue  may  be  stated  as  follows:  "Given 
laboratory  measurements  of  the  static  strength  of  rock,  what  are  the 
appropriate  values  to  be  used  in  simulating  the  response  of  rock  to 
an  underground  nuclear  explosion?"  The  uncertainties  associated 
with  answers  to  this  question  arise  because  laboratory  measurements 
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of  strength  vary  greatly  depending  whether  the  sample  is  wet  or  dry, 
fractured  or  intact,  and  on  the  strain  rate  at  which  the  test  is 
conducted. 

However,  there  is  a  minimum  strength  for  a  rock  mass;  i.e., 
the  limit  when  the  rock  is  completely  broken  into  blocks.  In  this 
case  Byerlee  (1979)  has  shown  that  the  strength  is  determined  by  the 
frictional  strength,  which  is 

x  *  0.85  a  (2.5) 

n 

where  t  is  the  shear  stress  and  a  the  normal  stress  applied  to 

n 

the  frictional  surface.  Equation  (2.5)  is  applicable  at  low  normal 
stresses  (0.2  Kb  <  o^  <  2.0  Kb).  At  higher  normal  stress,  the 
applicable  law  is 

T  -  0.5  +  0.6  On  (2.6) 

expressed  in  kilobars.  Equations  (2.5)  and  (2.6)  have  been  shown  to 
be  independent  of  lithology,  temperature  and  scale,  and  to  exhibit  a 
weak  rate  dependence  in  which  the  frictional  strength  decreases  a 
few  percent  per  decade  of  increase  in  sliding  velocity  (Scholz  and 
Engelder,  1976;  Dieterich,  1978). 

When  the  rock  contains  a  flu’d  within  its  pore  structure  at  an 
internal  pressure,  it  has  been  shown  (Garg  and  Nur,  1973)  that 

the  strength  of  the  rock  does  not  depend  only  on  the  externally 

applied  stress,  but  on  an  effective  stress  (o^  -  p^-).  For 

example,  when  fluid  is  present  in  the  rock.  Equations  (2.5)  and 

(2.o)  should  be  modifed  to  be 

t  »  0.85  (o^  -  pf)  (2.7) 

and 

x  «  0.5  +  0.6  (c^  -  pf)  (2.8) 

and  the  strength  is  said  to  obey  a  "law  of  effective  stress." 
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Even  the  synthesis  of  a  large  amount  of  data  as  represented  by 
Equations  (2.7)  and  (2.8)  involves  a  large  variation  in  strength,  t, 
depending  on  the  variation  of  pore  pressure,  p^,  with  the  stress 
state,  o.  One  might  assume  that  the  issue  could  be  resolved  by 
testing  rock  samples  obtained  from  ^ore  taken  at  the  site  of  the 
event.  In  general,  this  has  proven  to  be  an  unsatisfactory 
procedure  because  the  rock  sample  is  not  representative  of  the 
average  water  content  or  fracture  density  of  the  site,  and  because 
the  laboratory  induced  pore  pressure  does  not  correspond  to  that 
generated  in  situ  during  the  passage  of  the  shock  wave. 

Therefore,  due  to  the  wide  range  of  strength  values  available 
from  laboratory  tests,  the  only  way  to  decide  what  is  important  is 
to  try  to  match  ground  motion  data  from  explosive  sources,  compare 
the  strength  required  to  match  the  data  with  that  from  laboratory 
tests,  and  then  attempt  to  draw  general  conclusions  concerning  the 
types  of  laboratory  tests  most  appropriate  for  the  nuclear  test  rock 
environment.  That  is  the  approach  we  have  taken. 

The  conclusions  we  have  reached  regarding  the  strength 

required  to  simulate  explosive  induced  ground  motion  in  low  porosity 
brittle  rocks  (crystalline  igneous  rocks,  crystalline  metamorphic 
rocks,  and  well-consolidated  sedimentary  rocks)  and  high  porosity 
brittle  rocks  (poorly  consolidated  sedimentary  rocks  and  vesicular 
igneous  rocks)  are  as  follows: 

1.  The  strength  assumes  a  "high"  value  during  the  initial 

portion  of  the  stress  pulse  and  then  relaxes  to  a  "low" 
value. 

2.  The  "high"  value  appears  to  be  that  obtained  from 

laboratory  measurements  on  dry  rock  samples.  However, 
there  is  some  evidence  which  indicates  that  for  high 
strain  rates  (10^  -  10^/sec),  a  strain  rate  effect 
should  be  added  to  the  static  dry  strength  measurements. 

3.  The  low  value  appears  to  be  that  obtained  by  assuming 

that  an  effective  stress  law  is  operative  and  that  the 
pore  fluid  pressure,  p^,  is  equal  to  the  pressure  of 
the  externally  applied  stress,  P. 
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Apparently  the  first  to  recognize  the  importance  of  effective 
stress  for  matching  explosively  induced  ground  motion  data  were 
Cherry  and  Peterson  (1970).  They  assumed  that  the  wet  strength  was 
the  '‘equilibrium"  strength,  added  a  strain  rate  dependence  to  the 
dry  strength,  used  a  Maxwell  solid  formulation  to  relax  between  the 
dry  and  wet  strengths,  and  succeeded  in  matching  the  stress  history 
data  obtained  from  small  scale  experiments  in  blocks  of  grout. 
Figure  2.1  shows  the  differences  in  measured  strength  characteristic 
between  dry  and  wet  grout.  Figure  2.2  shows  the  comparison  between 
the  measured  and  calculated  stress  history,  6.5  centimeters  from  the 
charge,  in  grout  with  10  percent  air-void  porosity,  where  the 

calculation  was  performed  assuming  relaxation  between  the  strain 
rate  dependent  dry  strength  and  the  wet  strength. 

This  model  was  revised  by  Cherry  and  Rimer  (Bache,  et 
1975)  where  they  used  an  effective  stress  law  and  coupled  the 

Maxwell  solid  stress  relaxation  from  the  dry  to  the  saturated 
strength  with  the  removal  of  all  air-void  porosity  from  the  rock 

matrix. 

The  last  modification  to  the  strength  constitutive  model  has 
involved  replacing  the  Maxwell  solid  relaxation  with  a  direct 
calculation  of  the  change  in  pore  fluid  pressure  as  a  function  of 
the  parameter  a.  In  the  model,  we  assume  that  the  effective 
pressure,  P-p^,  varies  as  follows: 

P  -  Pf  ■  P  o  >  ag  (2.9a) 

p  “  pf  *  ~a ' ^  P  ae  >  a  >  1  (2.9b) 

P  -  Pf  -  0  a  -  1  (2.9c) 

where  a  and  P  are  defined  by  Equations  (2.1)  and  (2.3)  respectively, 
and 

a  is  the  value  of  a  at  P„. 
e  e 
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Figure  2. 


1.  Strength  of  DF-5A  grout  (Cherry  and  Peterson,  1970).  In  this 
figure  and  in  the  remainder  of  this  report,  Y  corresponds  to 
the  difference  between  the  maximum  and  minimum  principal  _ 
stresses  at  failure  under  triaxial  stress  conditions  and  P 
Is  half  the  sum  of  these  stresses.  For  a  definition  of  these 
parameters  In  terms  of  stress  invariants,  see  Cherry  and 
Peterson  (1970). 
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Figure  2.2.  Stress  history  in  DF-5A  grout  6.5  cm  from  high  explosive 
detonation  (Cherry  and  Peterson,  1970). 
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The  modifications  given  by  Equation  (2.9)  were  introduced  into 
the  model  both  to  improve  the  physics  in  the  model  so  that  effective 
stress  could  be  used  directly  in  the  calculation  of  strength,  and  to 
remove  the  strain  rate  dependence  of  the  Maxwell  solid  which 

obviates  the  cube  root  scaling  of  ground  motion  measurements  with 
explosion  yield.  The  material  strength  is  given  by 

Y  -  Y(F  -  pf)  (2.10) 

with  the  strength  function  determined  from  laboratory  tests  on  dry 
(pf*0)  rock  samples.  Equation  (2.10)  is  intended  to  be  a 

generalization  of  the  effective  stress  law  given  by  Equations  (2.7) 
and  (2.8).  The  reduction  in  strength  due  to  p^  is  introduced  into 
the  calculation  of  strength  via  the  effective  stress  constitutive 

model  given  by  Equations  (2.9)  and  (2.10). 

It  is  worth  emphasizing  that  the  effective  stress  constitutive 
model  makes  simple  and  straightforward  requirements  on  material 
properties  data,  with  the  critical  material  properties  being  the  dry 
strength,  Y(P),  and  the  pressure,  Pc,  at  which  all  air  voids  are 
removed.  Even  if  site  specific  measurments  of  these  mterial 
propc. w  es  are  not  available,  enough  laboratory  tests  have  been 
conducted  on  a  variety  of  rock  types  so  that  in  most  cases 
reasonable  estimates  of  these  properties  can  be  made,  e.g.,  Stephens 
and  Lilley  (1970),  Heard  (1970),  Byerlee  (1979). 

Our  first  application  of  this  model  involved  a  simulation  of 
ground  motion  data  recorded  from  the  60  KT  PILEDRI VER  event.  The 

material  properties  used  for  the  simulation  were  obtained  from  Heard 
(1970)  for  the  strength  of  dry,  fractured  granite,  and  from  Brace 
(1965)  for  the  pressure  (P  -  1  Kb)  at  which  all  void  space 
(fracture  porosity)  is  closed.  In  addition,  we  assumed  that  the 
site  contained  0.1  percent  air  void  porosity. 

Figure  2.3  and  2.4  compare  the  calculations  with  radial  ground 
motion  data  recorded  at  a  horizontal  range  of  204  meters.  Figures 
2.5  and  2.6  show  a  comparison  between  the  data  at  204  meters  and  a 
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Comparison  between  PILEDRIVER  displacement  data  (light  curve)  and  the 
effective  stress  law  calculation  (heavy  curve)  204  meters  from  the  source 
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calculation  in  which  the  effective  stress  law  was  not  used,  i.e. , 

the  pore  fluid  pressure,  p^,  in  Equation  (2.10)  was  assumed  to  be 

zero  for  all  strain  states  and,  therefore,  the  dry  granite  strength 
was  not  influenced  by  effective  stress.  The  difference  between 
these  two  calculations  is  quite  significant,  a  factor  of 
approximately  15  in  peak  displacement  for  instance,  with  the 
effective  stress  model  providing  good  agreement  with  the  data. 

In  addition,  the  final  cavity  radius  calculated  by  the 
effective  stress  model  was  40  meters,  also  a  good  agreement  with 
PILEDRIVER  cavity  measurements.  The  cavity  radius  from  the 
noneffective  stress  calculation  was  26  meters.  It  is  interesting 
that  the  difference  between  the  effective  stress  and  noneffective 
stress  cavity  radii  is  very  close  to  the  difference  between  NTS 
cavities  in  granite  and  those  reported  by  the  French  from  their 

granite  test  area  in  the  Sahara. 

The  effective  stress  model  was  used  by  Day,  Rimer  and  Cherry 
(1981)  to  perform  a  two-dimensional  ( ax i symmetric)  finite  difference 
simulation  of  the  PILEDRIVER  event  in  order  to  determine  the 
relation  between  spall  and  surface  wave  generation.  Figure  2.7 
shows  a  comparison  between  calculated  ground  motion  and  data 

recorded  from  the  event.  This  data  requires  a  two-dimensional 
simulation  because  of  the  influence  of  the  free  surface  on  the 
ground  motion. 

The  results  shown  in  Figures  2.3,  2.4  and  2.7  along  with  the 
comparison  with  cavity  radius  data  suggest  that  an  effective  stress 
law  is  able  to  provide  a  satisfactory  explanation  of  the  PILEDRIVER 
ground  motion  data.  Unfortunately,  this  type  of  data  from 
underground  nuclear  tests  is  not  frequently  obtained.  In  fact,  the 
PILEDRIVER  data  is  unique  in  terms  of  the  large  amount  of  data 
available  and  the  extent  of  both  surface  and  subsurface  coverage 
provided.  Therefore,  if  further  verifications  of  both  the 
air-filled  porosity  and  the  effective  stress  constitutive  models  are 
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Figure  2.7.  Vertical  velocity  time-histories  (positive  up)  obtained  from 
the  two-dimensional  simulation  of  Day,  Rimer  and  Cherry  (1981) 
compared  with  recorded  velocities  for  PILEDRIVER. 
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to  be  Obtained,  it  would  seem  logical  that  we  use  less  expensive, 
small  scale  explosive  tests.  Data  from  such  tests  has  recently 
become  available.  In  the  remainder  of  this  report,  we  show  the 
results  of  applying  these  models  to  this  data. 


III.  THE  GROUT  EXPERIMENTS 


As  part  of  the  Defense  Nuclear  Agency  (DNA)  late-time 
containment  program  for  underground  nuclear  testing,  SRI 
International  has,  since  1976,  been  conducting  small  scale 
laboratory  experiments  to  study  explosion  phenomenology  related  to 
the  containment  of  the  cavity  gas.  The  initial  experiments  involved 
casting  a  12  inch  diameter  sphere  of  either  2C4  or  low  density  LD2C4 
grout  around  a  lucite-encased  sphere  of  high  explosive  (3/8  gm  of 
PETN),  placing  this  grout  sphere  in  a  pressurized  water  tank  to 
simulate  overburden  pressure,  and  detonating  the  PETN.  While 
maintaining  overburden  pressure,  the  sphere  was  then  hydrofracture d 
from  the  explosively  formed  cavity  using  a  tube  emplaced  preshot. 
Hydrofracture  breakdown  pressures  from  these  tests  are  compared  with 
breakdown  pressures  from  precast  (unexploded)  spheres  to  obtain  an 
estimate  of  the  magnitude  of  the  explosively  formed  residual  stress 
fields.  Figure  3.1  shows  the  experimental  apparatus.  A  detailed 
description  of  these  experiments  and  a  discussion  of  the 
experimental  results  may  be  found  in  Cizek  and  Florence  (1981). 

Recently,  SRI  has  begun  to  emplace  particle  velocity  gauges 
inside  the  grout  spheres  (Cizek  and  Florence,  1981).  These  gauges 
consist  of  concentric  circular  current-carrying  loops  of  wire  cast 
symmetrically  about  the  charge.  A  magnetic  field  is  generated 
normal  to  the  plane  of  the  loops  by  passing  current  through  a  coil 
which  surrounds  the  sphere.  Charge  detonation  produces  radial 
motion  of  the  loops  that  cut  the  magnetic  flux  lines.  In  accordance 
with  Faraday's  law,  the  voltage  induced  in  each  conducting  loop  is 
proportional  to  the  particle  velocity.  Records  of  radial  particle 
velocity  versus  time  obtained  from  these  gauges  have  proved  to  be 
reproducible  from  shot  to  shot. 

In  addition  to  these  experiments,  DNA  has  supported  the 
measurement  of  material  properties  for  2C4  and  LD2C4  grout.  These 
measurements  have  been  made  at  SRI  and  Terra  Tek,  Inc.  (Cooley,  et 
al. ,  1982).  They  include  static  load-unload  tests  under  uniaxial 
strain  conditions,  hydrostatic  compression  tests,  triaxial 
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compression  tests  on  both  virgin  samples  and  on  samples  damaged  by 
hydrostatic  compression  to  4  Kbar  and  subsequent  decompression, 
ultrasonic  velocity  tests  and  porosity,  water  content  and  density 
measurements. 

Table  3.1  and  Figures  3.2,  3.3,  and  3.4  summarize  the  results 
of  these  tests.  The  data  given  in  the  table  and  in  the  figures  were 
used  in  the  effective  stress  simulation  of  the  experiments. 
However,  the  strength  data  of  Figure  3.4  were  modified  to  account 
for  the  fact  that  the  triaxial  tests  were  not  conducted  either  on 
dry  samples  or  in  an  environment  in  which  the  pore  fluid  pressure 
was  monitored  during  the  test. 

In  the  absence  of  strength  tests  on  dry  grout  samples,  we 
assumed  that  the  dry  failure  surface  was  simply  an  upward  extension 
of  the  sloped  portion  of  the  intact  failure  surface  shown  in  Figure 

3.4, 

Y  -  a  +  bV  .  (3.1) 

For  2C4  grout,  we  obtained  the  best  agreement  with  the  particle 
velocity  data  using 

a  -  0.09  Kb 

(3.2) 

b  *  1.6 

while  for  LD2C4  grout,  we  used 
a  -  0.055  Kb 

(3.3) 

b  *  l.o# 

Both  the  slopes  and  the  intercepts  given  by  Equations  (3.2)  and 

(3.3)  are  within  the  constraints  imposed  by  the  data  shown  in  Figure 

3.4.  The  grout  strength  is  now  modified  by  pore  fluid  pressure 
according  to 

Y  -  a  ♦  b(F  -  pf)  (3.4) 

where  p^  is  calculated  using  Equation  (2.9). 


TABLE  3.1 


MATERIAL  PROPERTIES  FOR  2C4  ROCK  MATCHING  GROUT 
AND  LD2C4  LOW  DENSITY  GROUT 


Property _ 2C4 _ LD2C4 


Density  (gm/cc)  2.205  1.921 

P-Wave  Velocity  (m/s)  3318.0  3128.0 

S-Wave  Velocity  (m/s)  1792.0  1629.0 

Air-Fi lied  Porosity  {%)  1.5  13.4 

Unconfined  Strength  (bar)  270  t  30  220  t  30 

Tensile  Strength  (bar)  36  ±10  31  +  10 
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surface  for  2C4  rock  matching  grout 


IV.  SIMULATION  OF  EXPERIMENTS  IN  2C4  AND  LD2C4  GROUT 


In  this  section,  we  present  three  simulations  of  high 

explosive  experiments  in  grout  and  compare  the  results  of  these 
simulations  with  particle  velocity  measurements  at  distances  of 
1.27  cm,  1.9  cm,  2.54  cm  and  4.0  cm  from  the  charge  center. 

The  first  simulation  used  the  material  properties  for  2C4 
grout  but  no  effective  stress  law,  and  the  strength  used  in  the 

simulation  was  that  denoted  in  Figure  3.4  as  the  "intact  failure 
surface."  Figures  4.1,  4.2,  4.3,  and  4.4  compare  the  particle 

velocities  calculated  from  the  simulation  with  those  measured  during 
the  experiment.  At  all  ranges,  the  calculated  velocity  pulse  is  too 
wide  during  the  positive  phase  and  too  narrow  during  the  negative 
phase. 

The  second  simulation  used  the  material  properties  for  2C4 
grout  and  included  the  modification  of  material  strength  due  to  pore 
fluid  pressure  (effective  stress)  given  by  Equation  (3.4).  The 

comparison  between  calculated  particle  velocities  and  those  measured 
are  shown  in  Figures  4.5,  4.6,  4.7  and  4.8.  At  all  ranges,  we 
obtain  excellent  agreement  with  measured  data  from  effective  stress 
law  simulation. 

The  consequences  of  these  two  simulations  (with  and  without 
effective  stress)  on  seismic  coupling  is  interesting  and  is  shown  by 
the  RVP  spectra  in  Figure  4.9.  The  RVP  spectrum  with  effective 
stress  is  peaked,  with  the  spectral  peak  about  a  factor  of  two 
larger  than  the  spectral  amplitude  at  low  frequencies.  The  RVP 
spectrum  without  effective  stress  is  flat  and  greater  than  the 
effective  stress  spectrum  at  low  frequencies  by  a  factor  of  1.5. 

Peaking  of  the  RVP  spectrum  is  a  direct  consequence  of  the 
width  and  amplitude  of  the  negative  portion  of  the  velocity  pulse, 
and  the  frequency  at  which  the  peak  occurs  is  approximately  equal  to 
the  period  of  the  pulse.  The  effective  stress  model  produces  the 
correct  shape  for  the  velocity  pulse  during  both  the  positive  and 
negative  phases  of  the  pulse  by  providing  high  material  strength 


Figure  4. 


(U>  Lo  io  an  4J  OjO  0£  7A  &0  8U0 

TIME  (SEC)  •10“ 


1.  Comparison  between  measured  velocities  in  2C4  grout 
and  those  calculated  from  a  noneffective  stress  law 
simulation  at  1.27  cm.  In  this  and  all  subsequent 
comparisons  the  solid  curves  are  calculated,  and  dashed 
curves  are  measured  velocities. 
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Figure  4.3.  Comparison  between  measured  velocities  in  2C4  grout 
and  noneffective  stress  law  simulation  at  2.54  cm. 
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Figure  4.4. 


Comparison  between  measured  velocities  in  2C4  grout 
and  noneffective  stress  law  simulation  at  4.0  an. 
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Figure 


4.6.  Comparisons  between  measured  velocities  in  2C4 
grout  and  effective  stress  law  simulation  at 
1 .9  cm. 
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Figure  4.7.  Comparisons  between  measured  velocities  in  2C4 
grout  and  effective  stress  law  simulation  at 
2.54  cm. 
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Figure  4.8.  Comparioons  between  measured  velocities  in  2C4 
grout  and  effective  stress  law  simulation  at 
4.0  cm. 
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3.  ROD  MODEL 


Figure  4.9.  Spectra  of  the  reduced  velocity  potentials  from  the 

noneffective  stress  and  the  effective  stress  simulations. 
Also  included  is  the  spectra  from  the  RDD  model  (Rimer 
and  Lee  (1982)) . 


34 


SCUBED 


during  the  initial  stages  of  pore  collapse  and  then  lowering  the 
strength  during  the  final  stages.  The  high  strength  retards  the 
velocities  during  the  positive  phase  while  the  low  strength  broadens 
the  negative  portion  of  the  pulse. 

The  third  simulation  used  the  material  properties  for  LD2C4 
grout  and  included  the  effective  stress  constitutive  model.  The 
comparison  between  the  calculated  particle  velocities  and  those 
measured  during  the  grout  experiments  are  shown  in  Figures  4.10, 
4.11,  4.12,  and  4.13.  The  agreement  between  the  calculated  results 
and  the  data  is  acceptable,  indicating  that  the  effective  stress 
model  is  applicable  to  a  wide  range  of  air  void  porosities. 
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gure  4.11.  Comparisons  between  measured  velocities  in  LD2C4 
grout  and  effective  stress  law  calculation  at 


Figure  4.12.  Comparisons  between  measured  velocities  in  LD2C4 
grout  and  effective  stress  law  calculation  at 
2.54  cm. 
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V.  UNIQUENESS 


Other  constitutive  models  can  be  constructed  which  provide  as 
good  (or  better)  a  fit  to  the  data  as  the  effective  stress  law 
results  shown  in  Section  IV.  Rimer  and  Lie  (1982)  have  developed 
such  a  model.  They  discuss  in  detail  a  complicated  rate  dependent 
and  damage  dependent  set  of  constitutive  models  (referred  to  as  the 
R00  model)  which  they  used  to  simulate  the  SRI  experiments  in  2C4 
and  LD2C4  grout.  The  R00  model  has  three  important  features:  a 
failure  surface  which  is  strain  rate  dependent  upon  loading;  a  rate 
dependent  stress  relaxation  from  this  failure  surface  to  a  damaged 
failure  surface  upon  unloading;  and  a  reduced  shear  modulus  upon 
loading.  Both  the  rate  of  deviatoric  stress  relaxation  and  the 
failure  strength  of  the  damaged  grout  are  made  functions  of  the 
amount  of  damage,  as  is  the  shear  modulus.  The  deviatoric  stress 
relaxation  to  the  damaged  failure  surface  is  accomplished  using  a 
Maxwell  solid  approach. 

This  model  successfully  matches  the  SRI  data  in  2C4  and  LD2C4 
grout  (Figures  5.1,  5.2,  5.3,  5.4,  5.5,  5.6,  5.7,  and  5.8),  and  the 
seismic  coupling  is  approximately  the  same  as  that  from  the 
effective  stress  law  (Figure  4.9).  However,  the  complicated  RDD 
model  requires  many  assumptions  regarding  material  properties  which 
are  certainly  not  easily  verified  (and  indeed  may  be  impossible  to 
verify)  from  laboratory  material  properties  tests. 


VELOCITY  (CM/SEC) 


Figure  5.1.  Comparison  between  measured  particle  velocities  from 
SRI  2C4  grout  test  273  and  results  of  RDD  model  cal¬ 
culation  at  a  radial  range  of  1.27  cm. 
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Flqure  5.2.  Comparisons  between  measured  velocities  from  2C4 
grout  tests  272  and  273  and  results  of  RDO  model 
calculation  at  1.9  cm. 


Figure  5.3.  Comparisons  between  measured  velocities  in  2C4 
grout  and  ROD  model  simulation  at  2.54  cm. 
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Figure  5. 


.  Comparisons  between  measured  velocities  in  2C4 
grout  and  ROD  model  simulation  at  4.0  cm. 
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Figure  5.5.  Comparisons  between  measured  particle  velocities 
from  SRI  L02C4  grout  tests  288  and  297  and  results 
of  ROD  model  calculation  at  a  radial  range  of 
1.27  cm. 
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Figure  5.6. 


Comparisons  between  measured  velocities  in  LD2C4 
grout  and  RDD  model  simulation  at  1.9  cm. 
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Figure 


5.7.  Comparisons  between  measured  velocities  in  LD2C4 
grout  and  RDD  model  calculation  at  2.54  cm. 
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Figure  5.8.  Comparisons  between  measured  velocities  in  LD2C4 
grout  and  ROD  model  calculation  at  4.0  cm. 
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VI.  CONCLUSIONS  ANO  RECOMMENDATIONS 


We  have  reviewed  the  current  status  of  two  specific 
constitutive  models  which  severely  affect  seismic  coupling.  These 
models  were  the  irreversible  collapse  of  air-filled  porosity  and  the 
modification  of  material  strength  due  to  pore  fluid  pressure,  i.e., 
effective  stress. 

We  have  shown  how  these  models  provide  a  very  simple, 
straightforward  explanation  of  ground  motion  data  over  a  wide  range 
of  yields  (61  KT,  nuclear  to  3/8  gram,  PETN)  and  air  void  porosities 
(0.1  percent  to  13.4  percent).  The  most  critical  material 
properties  required  by  these  models  for  near  field  ground  motion  and 
seismic  coupling  predictions  are  the  dry  strength,  the  pressure  at 
which  all  air  void  porosity  is  removed,  and  the  initial  air  void 
porosity. 

These  results  have  served  to  reenforce  our  confidence  in  the 
ability  of  these  models  to  predict  near  field  ground  motion  and 
seismic  coupling  effects.  We  recommend  intensive  continuation  of 
their  application,  not  only  to  explosion  sites  which  have  preshot 
measurements  of  material  properties  and/or  near  field  ground  motion 
measurements,  but  also  to  sites  where  explosion  yield  can  only  be 
estimated  from  far  field  seismic  measurements. 
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